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We review briefly how the thinking about the permeation of gases, especially CO2, 
across cell and artificial lipid membranes has evolved during the last 100 years. We then 
describe how the recent finding of a drastic effect of cholesterol on CO2 permeability 
of both biological and artificial membranes fundamentally alters the long-standing idea 
that CO2 — as well as other gases — permeates all membranes with great ease. This 
requires revision of the widely accepted paradigm that membranes never offer a serious 
diffusion resistance to CO2 or other gases. Earlier observations of "C02-impermeable 
membranes" can now be explained by the high cholesterol content of some membranes. 
Thus, cholesterol is a membrane component that nature can use to adapt membrane 
CO2 permeability to the functional needs of the cell. Since cholesterol serves many 
other cellular functions, it cannot be reduced indefinitely. We show, however, that cells 
that possess a high metabolic rate and/or a high rate of O2 and CO2 exchange, do 
require very high CO2 permeabilities that may not be achievable merely by reduction 
of membrane cholesterol. The article then discusses the alternative possibility of raising 
the CO2 permeability of a membrane by incorporating protein CO2 channels. The highly 
controversial issue of gas and CO2 channels is systematically and critically reviewed. It is 
concluded that a majority of the results considered to be reliable, is in favor of the concept 
of existence and functional relevance of protein gas channels. The effect of intracellular 
carbonic anhydrase, which has recently been proposed as an alternative mechanism to 
a membrane CO2 channel, is analysed quantitatively and the idea considered untenable. 
After a brief review of the knowledge on permeation of O2 and NO through membranes, 
we present a summary of the 18 0 method used to measure the CO2 permeability of 
membranes and discuss quantitatively critical questions that may be addressed to this 
method. 
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measurements 



THE CLASSICAL CONCEPT: LIPID BILAYERS AND CELL 
MEMBRANES ARE EXTREMELY PERMEABLE TO GASES 

It has been a long-standing assumption in biology that all cell 
membranes let gases pass extremely easily by virtue of the high 
gas solubility of lipophilic gases in the membranes' lipid phase. 
This view was formulated as a general principle of the perme- 
ation of substances across cell membranes as early as between 
1895 and 1901 by Ernest Overton (1895, 1896, 1899, 1901; see 
also review by Kleinzeller, 1997). Twenty years later, this same 
principle was shown to specifically hold for the gases CO2 and 
NH3 by Jacobs (1920a,b, 1922). The idea seemed to be supported 
by measurements of the gas permeabilities of artificial phospho- 
lipid membranes, which for O2, CO2, and N2 as well as other 
gases, turned out to be extremely high, i.e., considerably higher 
than the permeabilities of small uncharged polar molecules such 
as water, urea, and glycerol, not to speak of larger, uncharged 
polar molecules such as sugars or even of ions. This view has 
entered most textbooks as a well-established fact valid also for 
cell membranes (Alberts et al, 2002). In the case of CO2, which 



is highly lipophilic, but almost equally hydrophilic at the same 
time, the experimentally determined CO2 permeability values 
(PCO2) of lipid bilayers are around 1 cm/s (Gutknecht et al., 1977; 
Missner et al., 2008). How high this value is, may be appreciated 
by comparing it with the CI - — HCO^ permeability of the red 
cell membrane, which with a value of 1 • 10~ 4 -1 • 10~ 3 cm/s is 
3-4 orders of magnitude lower, although it is exceptionally high 
in comparison to other ion permeabilities. 

The concept of Overton has later been extended to the 
"solubility-diffusion mechanism" by Finkelstein (1976). With this 
modification, Finkelstein took into account that membrane per- 
meability may not only depend on the solubility in the membrane 
of the substance considered, but also on the diffusion resistance 
the substance encounters within the membrane. 

The transport rate of a substance across a membrane, dmldt, 
can be expressed by Fick's first diffusion law: 

dm/dt = — -D app - A ■ Ac w /d, or (1) 
dm/dt = —P-A- Ac w (2) 
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with D app representing the apparent membrane diffusion coef- 
ficient, A the diffusion area of the membrane, d the thickness of 
the membrane, P the membrane permeability, and Ac w the differ- 
ence of concentrations in the water phase immediately adjacent to 
two sides of the membrane. If the true diffusivity of the substance 
within the membrane, Dm, is to be considered, the concentration 
difference present within the membrane has to be expressed using 
the lipid/water partition coefficient Kp: 

dm/dt = —Dm ■ A ■ AcM/d = —Dm ■ A ■ Ac w ■ Kp/d, (3) 

and, from Equations (2) and (3): 

P = D M ■ Kp/d. (4) 

Thus, according to Finkelstein's concept, the membrane perme- 
ability will vary with the lipid-water partition coefficient, but will 
in addition depend on intramembrane diffusivity. 

While the intramembrane diffusion coefficients of gases have 
not been determined directly, the lipid-water partition coeffi- 
cients are fairly well known. This is shown for several gases of 
biological interest in Table 1. Firstly, it is apparent that the gases 
O2, NO, N2, and CO behave quite similar in terms of water 
and lipid solubilities as well as partition coefficients. CO2, on 
the other hand, differs greatly from all these gases with respect 
to solubility in water and in lipid. CO2 exhibits a 20-50 times 
higher solubility in water than the other four gases, and in oil 
an about 10-20 times higher solubility. This results in parti- 
tion coefficients that are only moderately smaller for CO2 than 
for the other gases: at most by a factor of 3-4 for the oil/water 
as well as for the phospholipid/water system. It is interesting 



Table 1 | Solubility behavior of gases in water and lipids. 




C0 2 


0 2 


NO 


N 2 


CO 


Solubility coefficient in H2O 
(25°C); [^imol/l/mmHg] 


45 a 


1.7 a 


2.6 b 


0.84 a 


1.3" 


Solubility coefficient in oil 
(25°C); [M-mol/l/mmHg] 


74 c 


6.5 C 




3.8 d 


5.3 e 


Oil-water partition coefficient 
(25°C) 


1.6 


3.8 




4.5 


4.1 


Phospholipid-water partition 
coefficient (25°C) 


0.95 f 


3. 98 


4.48 






Permeability in phospholipid 
membrane; [cm/s] 


~1 h 


4.1* 


4.6* 


2.8** 


2.6** 



"Battels et al. (1959); b Zacharia and Deen (2005); c Cottonseed oil, C0 2 
temperature-corrected, Lawrence et al. (1946); d Olive oil, Battino et al. (1984); 
e Olive oil at 3TC (temperature effect in oils is minor), Snedden et al. (1996); 
1 Egg lecithin liposomes, Simon and Gutknecht (1980); 9 Egg yolk phosphatidyl- 
choline liposomes, Moller et al. (2005); h Average from Gutknecht et al. (1977) 
and Missner et al. (2008). 'Predicted from Pco 2 phospholipid membrane and 
phospholipid-water partition coefficients. "Predicted from P C o 2 m phospholipid 
membrane and oil-water partition coefficients. 



to note that the gas partition coefficients in Table 1 in general 
appear quite similar for the oil/water and the phospholipid/water 
system. The last line of the table gives the measured Pco 2 m 
phospholipid membranes as mentioned above, together with 
the permeabilities predicted from Pco 2 by assuming that they 
vary proportionally to the respective partition coefficients. This 
implies the assumption that the intramembrane gas diffusion 
coefficients are identical for all five gases, which of course may not 
hold true. 

Under the condition of a constant Dm, the last line of Table 1 
suggests that all gas permeabilities should be rather similar, and 
all of them should be extremely high. It should be borne in mind, 
nevertheless, that CO2 is unique in Table 1 by possessing absolute 
solubilities in both water and lipid that are more than an order of 
magnitude greater than those of the other gases listed. 

A note of caution should be added concerning the above value 
of Pco 2 m lipid bilayers. In both determinations cited above, that 
of Gutknecht et al. (1977), who obtained Pco 2 = 0.35 cm/s, as 
well as that of Missner et al. (2008), who report Pqo 2 — 3-2 cm/s, 
it was attempted to eliminate the problem of unstirred layers 
by special precautions. As points of possible critique it should 
be mentioned that Gutknecht et al. (1977) showed that appar- 
ent bilayer CO2 permeability increased with rising concentrations 
[HCO3"] + [CO^] in the bathing solution (the latter causing 
increasing facilitation of CO2 diffusion in the unstirred layer), 
and finally reached a plateau value. However, this plateau does 
not necessarily prove that the diffusion resistance of the unstirred 
layer has completely vanished as was implied by the authors. 
Alternatively, the plateau might just indicate that facilitation of 
CO2 diffusion has reached a maximum, which could not be fur- 
ther augmented due to the increasing concentration of the two 
ions with their inhibitory effects on ionic mobility as well as 
on carbonic anhydrase activity. Thus, some doubt may remain 
how free of unstirred layer effects these measurements are. The 
technique employed by Missner et al. (2008) takes care of the 
unstirred layer problem by measuring pH gradients (and, derived 
from them, pco 2 gradients) within the unstirred layer in front 
of the bilayer. However, the extensive theory describing the pro- 
cess observed by these authors required fitting five parameters, 
one of them Pco 2 > to a single exponential curve of the type 
shown below (Figure 6A). Moreover, the fitting procedure did 
not seem to exhibit a clearly defined minimum of the sum of 
squared deviations, which prompted the authors to introduce a 
"penalty term" that prevented higher values of Pco 2 to emerge 
from the fitting procedure. Thus, both measurements leave some 
uncertainty about the correct value of Pco 2 m lipid bilayers. A 
further problem that is common to both approaches, results from 
the fact that in both studies planar lipid bilayers were employed, 
which were made from solutions of the lipids in n-decane. Bunce 
and Hider (1974) have shown that this may result in membranes 
whose final composition contains 30% decane. This is expected 
to lead to elevated Pco 2 values (a) because decane possesses a two 
times higher CO2 solubility (and partition coefficient) than phos- 
pholipid (Simon and Gutknecht, 1980), and (b) because decane 
may significantly increase membrane fluidity and thus Dm for 
C0 2 (Gutknecht et al., 1977). On the other hand, Gutknecht 
et al. (1977) used a 1:1 ratio of phospholipidxholesterol, which 
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is expected to reduce Pqo 2 ( see the following chapter). Missner 
et al. (2008), in contrast, used planar lipid bilayers with and with- 
out cholesterol, and observing no differences, they combined all 
results to calculate only an overall mean and SD. Nevertheless, 
in view of the results from other approaches to measure Pco 2 m 
phospholipid membranes, as discussed below, a value of ~1 cm/s 
as a crude estimate of Pco 2 m lipid bilayers appears reasonable. 

In contrast to the high CO2 permeabilities obtained by 
Gutknecht et al. (1977) and Missner et al. (2008), two other 
groups have obtained extremely low Pco 2 values from rapid reac- 
tion stopped-flow experiments with liposomes. Yang et al. (2000) 
measured by this technique the CO2 permeability of lipid vesi- 
cles without and with ~40% cholesterol, obtaining in both cases 
Pco 2 ~ 0.001 cm/s. Prasad et al. (1998) with the same technique 
also obtained a Pqo 2 of around 0.001 cm/s, no matter whether 
cholesterol was present or not. Yang et al. (2000) argued that the 
CO2 uptake process of vesicles in the stopped-flow apparatus is 
probably limited by unstirred layers. We may add that in some 
commercial stopped-flow apparatuses the mixing is less than per- 
fect, which would further increase the apparent unstirred layer. 
However, when cells or vesicles rather than just solutions are 
studied in a stopped-flow experiment, an additional complication 
arises: when the flow has stopped and cells or vesicles begin to take 
up some substance, then there is only diffusion, and no convec- 
tion, that can deliver the substance from the extracellular solution 
to the membrane. In other words, considerable gradients of this 
substance can develop in the extracellular/-vesicular space, and 
the standard assumption of a homogeneous concentration of this 
substance in the outside solution will not hold. This is not criti- 
cal, when the uptake process by cells or vesicles is slow, but in a 
case such as the extremely rapid uptake of CO2 into vesicles with 
a large CO2 binding capacity this can become a serious problem. 

Indeed, it can be shown that the process of CO2 uptake by 
a vesicle 150 nm in diameter with the buffer capacity of 25 mM 
HEPES takes <0.1 ms to reach 95% completion, when the chem- 
ical reaction inside the vesicles is infinitely fast and when there is 
no membrane diffusion resistance. 

This estimate is based on an approximation by Crank (1957) 
in combination with the use of the "effective solubility" within 
the vesicles, a concept applied to red cells by Thews (1960). 
Accordingly, effective solubility a' is the amount of CO2 plus 
HCO3" taken up by the intravesicular volume per increase in 
CO2 partial pressure. D in Crank's treatment is then replaced 
by D' = D ■ a/a! , where a is the classical solubility of molecular 
CO2. Crank's expression for 95% completion of the CO2 uptake 
by a sphere then becomes: 

D' ■ t/r 2 = 0.6, 

where t is the time required to reach 95% completion, and r is the 
radius of the sphere. For 25°C, with D = 1.8 • 10~ 5 cm 2 /s, a = 
9.8 • 0~ 4 cm 3 /cm 3 /mmHg, a' = 53 • 10~ 4 cm 3 /cm 3 /mmHg (this 
number being calculated for apco 2 jump from 1 to 40 mmHg), t 
from the above equation turns out to be 1 • 10~ 5 s = 0.01 ms. 

Thus, without a drastic diffusion resistance of the vesicle mem- 
brane to CO2, a stopped-flow apparatus with a dead time of a few 
milliseconds, and perhaps further problems limiting the speed of 



the process, cannot be expected to capture the kinetics of CO2 
uptake, provided that substantial intravesicular carbonic anhy- 
drase activity is present and the CO2 hydration kinetics assumes a 
half-time of < 1 ms, as was the case in both studies. We conclude 
therefore that the low Pqo 2 values reported by Prasad etal. (1998) 
and Yang et al. (2000) are likely to be essentially determined by 
technical problems of the stopped-flow experiment. 

The third approach that has been applied to determine Pqo 2 of 
artificial liposomes is the mass spectrometric ls O exchange tech- 
nique, which measures very fast processes on an extremely slow 
time scale (Itel et al, 2012, and see last section of this article). 
Unstirred layers have been shown to play a minor role (Endeward 
and Gros, 2009). In the absence of cholesterol, this method 
gives, as described below, a Pco 2 > 0.16 cm/s, where this num- 
ber reflects the upper limit of detectability of the technique for 
vesicles. Summarizing the discussion of this paragraph, it appears 
to be the best guess at present that the true Pco 2 of cholesterol- 
free lipid bilayers is certainly >0.16 and probably around 1 cm/s. 
Figure 1 gives a summary of the Pco 2 values obtained for arti- 
ficial lipid bilayers, by the stopped-flow technique using vesicles, 
by observing steady state CO2 transfer across planar lipid bilay- 
ers (Gutknecht et al, 1977; Missner et al, 2008), or by 18 0 mass 
spectrometry applied to vesicles (Itel et al., 2012). The drastic dif- 
ference between results from stopped-flow experiments and those 
from the other techniques is obvious. 

Recent work has called into question the classical simple view 
that cellular gas permeability is extremely high and reflects essen- 
tially the lipid/water partition behavior of the gas, a view that has 
recently, for example, been re-stated by Missner and Pohl (2009). 
In the case of CO2 several alternative mechanisms are currently 
being discussed. In the following sections, we will critically discuss 
the concepts presently proposed. 



Stopped-flow 



3 0 



1 -; 



E 0.1 
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o 0.01 
0_ 



0.001 



0.0001 



Artificial Liposomes 



i 4 p pH micro- 
electrode 



Black Lipid 
Membranes 



FIGURE 1 I Overview of C0 2 permeabilities determined for artificial 
liposomes or planar lipid bilayers using various techniques. References 
see text. 18 0 refers to the 18 0 exchange technique explained in the last 
section of this article. 14 C-labeled C0 2 was used by Gutknecht et al. (1977) 
to measure steady state transfer of C0 2 . Missner et al. (2008) have used a 
pH microelectrode to measure the gradient of pH, dpH/dx, in the unstirred 
layer in front of the membrane to be studied. From this gradient the flux of 
C0 2 across the diffusion layer was derived as explained in more detail in 
section V of this article. 
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EFFECT OF CHOLESTEROL ON C0 2 PERMEABILITY 
OF MEMBRANES 

Missner et al. (2008) as well as Prasad et al. (1998) and Yang et al. 
(2000) have observed no effect of cholesterol on Pqo 2 °f artificial 
lipid membranes. Possible limitations of Missner's experimental 
approach, e.g., the presence of decane in their planar bilayers and 
problems with their fitting procedure, have been discussed above. 
Likewise, it appears probable that an indequate time resolution of 
the stopped-flow experiments employed by Prasad et al. (1998) 
and by Yang et al. (2000) slowed down the observed kinetics so 
drastically that any moderate changes in true membrane perme- 
ability could not become apparent. In view of this it is interesting 
that Itel et al. (2012) have recently reported that cholesterol has a 
most drastic effect on the CO2 permeability of phospholipid vesi- 
cles. This is illustrated in Figure 2, which is reproduced from their 
paper. It can be seen that between 0 and 30 mol% cholesterol con- 
tent of the vesicles Pco 2 falls from >0.16 cm/s to ~0.02 cm/s, i.e., 
by at least one order of magnitude. If Pco 2 at 0 0/0 cholesterol were 
1 cm/s as given in Table 1, this fall would even occur over two 
orders of magnitude. Between 30 and 70 mol% cholesterol then, 
Pqo 2 decreases by another order of magnitude to ~0.0024 cm/s. 
It was concluded that cholesterol between 0 and 70% can reduce 
vesicle Pco 2 by 2-3 orders of magnitude. This drastic effect of 
cholesterol does not seem to be due to a decrease in CO2 solubil- 
ity in the lipid phase, as Simon and Gutknecht (1980) report only 
a 25% reduction of CO2 solubility when 50% cholesterol is added 
to egg lecithin in liposomes. Thus, cholesterol is likely to cause a 
significant reduction of intramembrane CO2 diffusivity Dm [see 
Equation (4)]. 

It is very important to note that the results of Figure 2 are 
in excellent agreement with the molecular dynamics simula- 
tions of Hub et al. (2010), who also predict that cholesterol can 



reduce the CO2 permeability of lipid bilayers by several orders 
of magnitude. Itel et al. (2012) have confirmed the quantita- 
tive findings of Figure 2 by a qualitative approach applied to 
intact cells. Decreasing the cholesterol content of MDCK cells by 
cholesterol depletion with |3-cyclodextrin caused a dramatic ^40- 
fold increase in CO2 permeability. Inversely, enrichment of these 
membranes by cholesterol-preloaded P-cyclodextrin reduced 
Pco 2 to 1/3 of normal MDCK cells. 

Does the finding shown in Figure 2 for artificial vesicles also 
apply to cell membranes? Itel et al. (2012) have shown that the 
CO2 permeability of cell membranes that are devoid of protein 
gas channels, as will be discussed below in chapter VI of this 
article, can be fairly precisely predicted from their cholesterol con- 
tent, when the dotted regression line of Figure 2 is used. Thus, 
it was concluded that also in cell membranes, in the absence of 
gas channels, cholesterol is the decisive determinant of CO2 per- 
meability and therefore is much more important than the exact 
phospholipid composition and membrane protein content. Since 
the membranes of many cells possess a cholesterol content of 
around 40%, a Pqo 2 of close to 0.01 cm/s is predicted from the 
regression line of Figure 2. Indeed, Itel et al. (2012) report for all 
of four cell types, namely two gas channel-free cell lines (MDCK 
and tsA201 cells), the red cell membrane without functional gas 
channels, and the presumably gas channel-free basolateral mem- 
brane of proximal colon epithelial cells, Pco 2 values of around 
0.01 cm/s. In conclusion, the permeability of many cell mem- 
branes is two orders of magnitude lower than the Pco 2 value of 
cholesterol-free artificial lipid membranes, which has long been 
believed to be representative for cell membranes too. 

"CO2-IMPERMEABLE" MEMBRANES 

A very early report of a cell membrane poorly permeable to 
a gas describes the observation of Kikeri et al. (1989) that the 
apical membranes in the medullary thick ascending limb of 
Henle in mouse kidney are "impermeable" to NH3. A few years 
later several instances of membranes impermeable to CO2 were 
reported, always apical membranes of epithelial cells. Waisbren 
et al. (1994) observed an apparent impermeability for CO2 of 
the apical epithelial chief and parietal cell membranes in iso- 
lated perfused gastric glands. Hasselblatt et al. (2000) reported the 
same property for the apical epithelial membranes of in vitro per- 
fused rat colonic crypts. Endeward and Gros (2005) then made 
an analogous observation on the surface cells of guinea pig colon 
epithelium, and for the first time they were able to quantitate this 
"impermeability"; they determined by ls O exchange mass spec- 
trometry the CO2 permeability of the apical proximal as well as 
distal colonic epithelial membrane to be ~0.001 cm/s. This was 
10 times lower than the CO2 permeability of several gas channel- 
free cell membranes as discussed above. The mechanism of the 
very low CO2 (and NH3) permeability of several apical epithelial 
membranes was entirely unclear, and only recently the results of 
Itel et al. (2012) have elucidated the cause of this property. Quite 
fortunately, Meyer zu Diittingdorf et al. (1999) have been suc- 
cessful in obtaining a good separation of apical and basolateral 
membranes of guinea pig colon epithelium and they performed 
an analysis of their lipid composition. They found a cholesterol 
content of 42% in the basolateral membrane of the proximal 



>0.16 cm/s 
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FIGURE 2 | Dependence of the C0 2 permeability of phospholipid 
vesicles on their cholesterol content. The latter is given on the x-axis in 
mol% per total membrane lipid. Vesicle phospholipid composition was 
phospatidylcholine:phosphatidylserine = 8:2. P cc . 2 was determined by the 
mass spectrometric 18 0 exchange technique. Dotted line is the linear 
regression line for all data points between 30 and 70% cholesterol. Pco 2 
values between 0 and 17% cholesterol were above the upper limit of 
detectability of the method of 0.16cm/s. Figure reproduced from Itel et al. 
(2012), with permission. 
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FIGURE 3 | Scheme of the model used to estimate the effect of 
membrane C0 2 permeability on the release of C0 2 from a 
cardiomyocyte. CO2 diffuses from the CO2 source located in the center of 
the cell first through the cytoplasm, and then through the plasma 
membrane, into the blood capillary. 



colon, which is what many cells have, but as much as 77% choles- 
terol in the apical membrane. Using the regression line of Figure 2 
one predicts for this cholesterol content a CO2 permeability close 
to 0.001 cm/s. This means that the very low permeability of this 
membrane can be satisfactorily explained simply by its cholesterol 
content. We conclude that cells can regulate their CO2 permeabil- 
ity over a wide range (and likely their NH3 permeability as well) 
by altering the cholesterol content of their membrane. 

The physiological advantage of making cell membranes imper- 
meable to gas is obvious in the case of NH3, where the potentially 
very high partial pressures in the lumen of the stomach as well 
as the colon would be toxic to the epithelial cells if the NH3 had 
access to their interior. Likewise, the potentially very high CO2 
partial pressures (up to 0.5 atm) in both lumina would constitute 
a severe acid load for the cells, if the CO2 were able to enter them 
(Endeward and Gros, 2005). 

WHICH C0 2 PERMEABILITY IS NECESSARY FOR RESPIRING 
CELLS? 

Itel et al. (2012) have reported Pco 2 values of 0.01 cm/s for 
two cell lines. Cells in culture usually exhibit a low oxidative 
metabolism, as they are not in an active state and often cover part 
of their energy requirement by anaerobic metabolism. Would a 
permeability of 0.01 cm/s suffice, for example, for highly active 
cells such as cardiomyocytes in vivo under conditions of heavy 
exercise? We present in the following calculations an answer to 
this question, employing a simplified model of CO2 release from 
a cardiomyocyte. The cardiomyocyte is considered as a cuboid, 
which in conjunction with further simplification allows us to 
model CO2 release on the basis of diffusion across a plane sheet. 
As the radius of human cardiomyocytes is 7 \im, we approximate 
the geometry of the cell by a plane sheet of 14 [im thickness, and 
perform the calculation for the half-thickness of the cuboid of 
7 |xm (Armstrong et al., 1998; Endeward et al., 2010; Endeward, 
2012). In Figure 3, the upper border of the cytoplasmic space then 
represents the center of the cell, where for simplicity the source of 
CO2, i.e., the mitochondria, is assumed to be located. From there 
the CO2 diffuses through the cytoplasm over 7 u,m, and then dif- 
fuses through the cell membrane before it is taken up by the blood 
in the capillary. We consider steady state conditions, so the CO2 
partial pressures in the center of the cell and in the capillary are 
constant. 

What is a realistic pco 2 gradient within the cell? It is gen- 
erally accepted that pco 2 gradients in tissues — in contrast to 
po 2 gradients — are extremely small, but no information on their 
precise value is available. Therefore, we use here the following 
approach. For the condition of heavy exercise, the heart can 
assume an O2 consumption, and about identical CO2 produc- 
tion, of 0.4 ml gas/g tissue/min (Endeward et al., 2010). With the 
known cytoplasmic CO2 diffusion coefficient, 1.3 • 10~ 5 cm 2 /s, 
and the diffusion path given in Figure 3, we calculate the par- 
tial pressure difference, Apco 2 > necessary to produce a flow of 
CO2 of 0.4 ml CGVg tissue/min from CO2 source to CO2 sink. 
In this calculation, the cell membrane is assumed to offer no dif- 
fusion resistance to CO2 at all. This yields a necessary Apco 2 of 
0.35 mmHg as indicated in Figure 3. This is a reasonable number 
in view of the expected small CO2 gradients in tissues. We then 



add the diffusion resistance of the membrane to that of the 7 \im 
thick cytoplasmic layer, where the latter amounts to 54s/cm. 
Using Fick's first law of diffusion, the rate of CO2 production, 
Vco 2 i and of diffusion from the CO2 source to the capillary, is 
then given by: 

V C02 = (d/D C02 + l/Pcc^r 1 • A • a C o 2 • Ap C o 2 /(d • A), (5) 

where the terms in the first brackets represent the sum of the resis- 
tances of diffusion through the cytoplasm and across the mem- 
brane, d is the cytoplasmic diffusion path, Dco 2 the cytoplasmic 
CO2 diffusion coefficient, Pco 2 membrane CO2 permeability, A 
the unit diffusion area (1 cm 2 ), aco 2 the solubility of CO2 in tis- 
sue (7.2 • 10~ 4 mm Hg _1 ), and the product d-A is the volume 
of cells in which Vco 2 is produced. This equation leads to the 
curve shown in Figure 4. At Pco 2 = 00 the value of Vco 2 is, as 
explained above, the maximal rate of cardiac CO2 elimination 
of 0.4ml/min/g. At constant Apco 2 > the CO2 production rate 
decreases with decreasing Pco 2 ■ While a Pqo 2 of 1 cm/s has only a 
small limiting effect on the rate of CO2 release, the Pco 2 found in 
cells in culture, 0.01 cm/s, would reduce CO2 production/release 
to 1/3 of the value in the absence of a membrane resistance. 
Pqo 2 = 0.001 cm/s would even reduce it to 1/20. Figure 4 shows 
that Pco 2 values <0.1 cm/s become increasingly critical for cellu- 
lar CO2 release, an effect that can only be overcome by raising 
intracellular pco 2 > and thus Apco 2 [see Equation (5)]. It can 
also be deduced from Equation (5) that a lower rate of CO2 
production, for example that of 0.003 ml/g/min as it occurs in 
resting skeletal muscle, which is only 1% of the value in maxi- 
mally working heart, can at identical or lower Apco 2 easily be 
achieved without any limitation with a Pco 2 of 0.01 cm/s, even 
with_Pco 2 = 0.001 cm/s a Apco 2 of 1/5 of 0.35 mmHg would suf- 
fice. Thus, it is clear that, while Pqo 2 = 0-01 cm/s is more than 
adequate for cells with low oxidative metabolism, a >10 times 
higher Pco 2 would seem important for cells with an extremely 
high oxygen consumption. It remains to be tested, whether such 
cells indeed exhibit CO2 permeabilities of that size. 

A different approach to the same question has been reported 
by Endeward et al. (2008), who evaluated the effect of membrane 
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FIGURE 4 | Influence of membrane C0 2 permeability on C0 2 release 
from a maximally metabolizing cardiomyocyte at a given partial 
pressure gradient of 0.35 mmHg over the half-thickness of the cell. 



CO2 permeability on the kinetics of gas exchange of red 
cells. Human red cells normally possess a CO2 permeability of 
0.15 cm/s. With this value, the time required for 95% completion 
of the release of CO2 as it occurs in the lung was calculated to 
be 110 ms, which may be compared to the capillary transit time 
in the lung of 700 ms. If red cells without functional gas chan- 
nels, which have a Pco 2 °f 0-01 cm/s (Endeward et al, 2008), are 
considered, the required time rises to 1000 ms, thus exceeding the 
available transit time in the lung. The limitation of CO2 release 
is expected to be more pronounced under conditions of exercise, 
when the transit time can fall to 350 ms. The problem caused by 
such a low _Pco 2 can on ly be overcome by increasing the pco 2 gra- 
dient across the alveolar diffusion barrier, and it is clear that the 
normal high CO2 permeability of red cells will ensure a maximally 
economical CO2 release in the lung. 

We conclude from these two approaches that optimal CO2 
release in cells with very high rates of gas exchange, like cardiomy- 
ocytes and red cells, requires high CO2 permeabilities >0.1 cm/s. 

IS INTRACELLULAR CARBONIC ANHYDRASE RESPONSIBLE 
FOR THE HIGH C0 2 PERMEABILITY OF SOME CELLS? 

Recently, Zocher et al. (2012) proposed that the low perme- 
ability of the urinary bladder with a Pco 2 of ~0.00016cm/s is 
not due to the uroplakin of the bladder epithelium but to the 
absence of carbonic anhydrase in that epithelium, a hypothesis 
constituting a basically new concept of membrane CO2 trans- 
fer. Measurements of CO2 transfer across monolayers of MDCK 
cells, which are known to possess intracellular carbonic anhy- 
drase, seemed to confirm this hypothesis: (1) when Zocher et al. 
(2012) added the carbonic anhydrase inhibitor acetazolamide to 
the solutions to inhibit all carbonic anhydrase (intra- as well as 
extracellular), the CO2 flux across the epithelium was very low, 
(2) it increased moderately when no acetazolamide was added, 



and (3) it increased 14-fold in comparison to experiment (1) 
when they added a very high activity of carbonic anhydrase (as 
produced by 2 mg CA/ml) to the solutions of the bathing medium 
of the epithelium. They concluded from this that "carbonic anhy- 
drase expression may uniquely regulate the tightness of biological 
membranes to CO2" and that their results "rule out that aqua- 
porins significantly contribute to the maintenance of acid base 
homoeostasis." Thus, these authors propose a biological mecha- 
nism that would raise the CO2 permeability of a membrane not by 
reducing membrane cholesterol, as discussed above, nor by inser- 
tion protein gas channels into the membrane, as we will discuss 
below. 

In order to understand the mechanism underlying the results 
of Zocher et al. (2012), we modeled the experimental situation 
used in their experiments with MDCK monolayers. The princi- 
ple of their measurements was to expose a layer of epithelium to 
a pco 2 gradient, and then measure the pH gradient that exists 
after steady-state has been reached in the unstirred layer between 
the front side of the epithelium (see Figure 5) and the bulk solu- 
tion in the downstream compartment. They do this by moving 
a pH microelectrode from the bulk solution towards the surface 
of the epithelium and recording pH as a function of the distance 
from the epithelial surface. Figure 6A shows examples in a figure 
reproduced from their paper (Zocher et al., 2012). Up to this step 
their procedure is identical to that used by Missner et al. (2008), 
which has been discussed above. However, they avoid the exten- 
sive theory required by Missner's approach by using an elegant 
approximation in the evaluation of these experiments. Zocher 
et al. (2012) use the first linear parts of the curves of Figure 6A to 
calculate the flux of CO2 across the epithelium, roughly approxi- 
mated by the HCO^ flux = (buffer-facilitated) H + flux through 
the unstirred layer in front of the epithelium: 

transepithelial CO2 flux ~ buffer-facilitated proton flux in 

unstirred layer 

= D B • BF • ApH/Ax, (6) 

where Dg is the diffusion coefficient of the buffer present in the 
extracellular solution, BF is the buffer factor provided by this 
buffer (in mM/ApH), and ApH/Ax is the linear initial slope of 
the curves in Figure 6A. 

Figure 5 illustrates the experimental situation employed by 
Zocher et al. (2012) as well as the principle of the present math- 
ematical model. Our model considers two unstirred layers on 
either side of the epithelial layer, plus the intracellular compart- 
ment. At the left hand side of the back-side unstirred layer and 
at the right-hand side of the front-side unstirred layer, the con- 
ditions of the bulk extracellular solution are fixed: pH = 7.4 in 
NaCl-HEPES buffer with 25 mM bicarbonate. pco 2 a t the back 
side was then increased abruptly in the model by raising [HCO^ ] 
from 25 to 50 mM in the left-hand side bulk solution, and sub- 
sequently pH at the front-side became acidified by the ensuing 
flux of CO2 across the epithelium. In all three compartments 
we modeled the transport of CO2, HCO^, H + , and HEPES 
(as a H + -carrier) by considering all diffusion and reaction pro- 
cesses involved. The membranes were considered impermeable 
to HCO3" and were assigned a CO2 permeability of 0.01 cm/s 
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FIGURE 5 | Scheme of the mathematical model used to simulate the 
CO2 flux measurements of Zocher et al. (2012). MDCK cell monolayers 
were exposed to an extracellular solution on the apical (back) side as well as 



on the basal (front) side. The flow is determined by measuring the pH 
gradient with a pH microelectrode moving across the front-side unstirred 
layer and by inserting this pH gradient into Equation (6). 




FIGURE 6 | pH gradients in the unstirred layer in front of the epithelial 
layer of MDCK cells. Acidification relative to the pH in the bulk 
extracellular solution occurs due to the CO2 influx from the epithelium. (A) 



Measured pH profiles reproduced with modification from Figure 7A of 
Zocher et al. (2012), with CC-BY license. Acet. = acetazolamide, (B) 
calculated pH profiles using the model given in Figure 5. 



(Itel et al, 2012). The differential equations were solved by a finite 
difference method as implemented in MATLAB 2008b. Each of 
the three compartments of Figure 5 was divided into 50 volume 
elements for this calculation. 

Figure 6A shows the experimental results obtained by Zocher 
et al. (2012). As described above, these authors used the initial 
slopes of the curves in Figure 6A to calculate the approximate 
CO2 fluxes occurring across the monolayer of MDCK cells. From 
these slopes as read from the curves of Figure 6A, the flux result- 
ing in curve (3; red) is estimated to be 6.3 times greater than 
the flux leading to curve (1; blue). This was interpreted to mean 
that in the presence of intra- (as well as extracellular) CA the 
epithelium possesses a higher CO2 permeability, and when no 
extracellular CA is added and all the endogenous CA of the 



MDCK cells is inhibited by acetazolamide, the CO2 permeability 
is 6.3 times lower. In the absence of the CA inhibitor, the curve (2; 
purple) assumes an intermediate position due to the endogenous 
CA activity of the cells. The conclusion was that carbonic anhy- 
drase, intracellular and ostensibly also extracellular, substantially 
increases CO2 permeability of the cell layer, which was attributed 
to an increase in cell membrane Pco 2 by Zocher et al. (2012). 

Figure 6B shows the results of the present model calcula- 
tions. In contrast to Figure 6A, we consider the unstirred layer 
only, without the bulk solution, and assume the unstirred layer 
to extend over 130 u,m in front of the MDCK cell layer. At the 
outer end of this unstirred layer (right hand side in Figure 6B), 
we hold the conditions of the extracellular bulk solution fixed 
(pH = 7.4, [HCO3"] = 25 mM). Reassuringly, our model predicts 
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a similar about eight times greater initial slope of curve (3; red) 
than of curve (1; blue), where in analogy to Figure 6A, curve (3) 
was calculated for an extracellular CA activity Aca>£ of 10,000 
and an intracellular activity Aca>1 of 800, which is realistic for 
MDCK cells (Al-Samir et al., 2013). In contrast, curve (1) was 
calculated for the complete absence of CA in the extra- as well 
as the intracellular space, which corresponds with the addition 
of 1 mM acetazolamide in Zocher's experiments. We conclude, 
that the model of Figure 5 nicely explains curves (1) and (3) of 
their experiment as it shows a similar difference in pH slopes 
to that observed by them. What is the mechanism of the inter- 
mediate (~4-fold) increase in slope of their curve (2) over that 
of curve (1)? Curve (2) of Figure 6A was obtained without any 
additions at all. A four-fold acceleration of the intracellular CO2 
flux by facilitated CO2 diffusion due to Aca>1 is unlikely in view 
of the findings that facilitated diffusion inside red cells as well 
as inside skeletal muscle (Gros and Moll, 1972; Kawashiro and 
Scheid, 1976) accelerates intracellular CO2 transport by no more 
than 60-80%. We have speculated therefore, that the slope of 
curve (2) in Figure 6A may be increased due to the combined 
action of intracellular and a small amount of membrane-bound 
extracellular CA, which both occur in MDCK cells (Pfaller et al., 
1989). Modeling this with an intracellular CA activity Aca.i of 
800 and an assumed membrane-bound surface-associated extra- 
cellular activity Aca>s of 100, we obtain curve (2) in Figure 6B, 
whose slope is indeed intermediate between those of curves (1) 
and (3). We note that in this model the activity Aca,s = 100 
is only applied in the one volume element immediately adja- 
cent to the cell membrane, while CA activity is assumed to be 
absent from all other volume elements of the unstirred layers. 
In conclusion, our model satisfactorily explains all experimen- 
tal results of Zocher et al. (2012) as shown in Figure 6A for 
MDCK cells. 

Nevertheless, the present model does not confirm the con- 
clusion of Zocher et al. (2012) that the flux of CO2 across the 
cell layer depends strongly on CA activity. This is due to the 
basal assumption of their method, which is that the CO2 flux 
leaving the cell layer and entering the right-hand side extracellular 
unstirred layer compartment (Figure 5), can in all cases be 
deduced from the surface pH gradient and the proton flux derived 



from it by Equation (6). Such a calculation is only a reasonable 
approximation, however, if the extracellular carbonic anhydrase 
activity is very high and therefore the CO2-HCO3 -H + -system is 
in near-complete chemical equilibrium. A CAII (the isoform II 
of CA) concentration of 2 mg/ml in the NaCl-HEPES solution, 
which was used experimentally, is more than the CAII concen- 
tration present in human red cells, and therefore the activity 
of 10,000 as used in our model is a minimum estimate of the 
CA activity in that solution (Wistrand, 1981; Endeward et al., 
2006a,b; Al-Samir et al., 2013). Thus, the assumption of near- 
perfect chemical equilibrium is certainly justified for curves (3) 
in Figures 6A,B- The CO2 fluxes derived from these slopes should 
be correct (and in Figure 6B indeed come close to the true ones, 
as derived from the model and given in Table 2, for the case 
of curve 3). However, in the absence of any CA activity, as in 
curve ( 1 ) of Figure 6A, and also in the presence of only very lim- 
ited extracellular CA activity as is likely the case in curve (2) of 
Figure 6A, it is not permissible to derive CO2 fluxes from the 
pH gradient. When the speed of the CO2 hydration reaction is 
slow, the CO2 leaving the cell layer will not react fast enough to 
produce the amount of protons that would be produced in the 
presence of high CA activity, and the C02-HCO^-H + -system 
will remain far from chemical equilibrium. Thus, deriving CO2 
fluxes from the extracellular pH profile using Equation (6) in 
the manner described, will lead to a substantial underestima- 
tion of these fluxes in the cases of curve (2), and even more 
so of curve (1). This is illustrated in Table 2, which gives the 
CO2 fluxes calculated from the mathematical model of Figure 5. 
It is apparent that that these fluxes do not vary by a factor of 
8, as the initial pH slopes do, or even 14, which is the aver- 
age factor of increase in CO2 flux given by Zocher et al. (2012). 
Instead, while the fluxes do show some minor positive correla- 
tion with extracellular CA activity, this effect does not exceed an 
increase in CO2 flux by more than 25%. This small effect is due 
to the acceleration of CO2 hydration on the membrane surface, 
which will remove the CO2 just having crossed the membrane 
more rapidly from the extracellular surface region of the cell 
membrane, thus maintaining a somewhat larger gradient of CO2 
across the membrane. This of course does not reflect any alter- 
ation of the CO2 permeability properties of cellular membranes. 



Table 2 | Summary of the conditions in the experiments with MDCK monolayers by Zocher et al. (2012) and in the model of Figure 5. 





Experimental 
condition 


CA activities used in the model, and 
likely present also in experiment 


Initial pH slope in front-side unstirred 
layer (from model. Figure 6B) (pH/cm) 


True CO2 flux derived from model 
(Figure 6B) (mmoles/cm 2 /s) 






AcaJ 


AcA.e A CA ,s 






Curve 1 


1 mM 






1.6 


2.0- 10" 6 




acetazolamide 










Curve 2 


No additions 


800 


100 


5.1 


2.2- 10" 6 


Curve 3 


2 mg/ml CA in 


800 


10,000 


12.6 


2.5- 10" 6 




extracelluar buffer 











A C aJ intracellular CA activity, A CA ,e extracellular activity in the bathing solutions, A C a,s assumed value of membrane-bound, cell surface-associated CA activity, which 
is active only in a small volume element immediately adjacent to the cell membrane. The CO2 fluxes in the rightmost column are derived from the mathematical 
model, and for curves 7 and 2 deviate substantially from those obtained from the slopes of the curves in Figure SB (2nd column from right). If it were permissible 
to consider the slopes of the 2nd column from right as proportional to the CO2 fluxes, an 8-fold decrease in flux under the condition of curve 1 would be predicted 
in comparison to the condition of curve 3. 
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We conclude that carbonic anhydrase is not involved in regulat- 
ing the CO2 permeability of cell membranes per se, and other 
mechanisms must be responsible for the considerable variations 
of CO2 permeability found between different membranes (Itel 
etal., 2012). 

PROTEIN CHANNELS FOR MOLECULAR C0 2 IN 
MEMBRANES OF HIGH C0 2 PERMEABILITY 

While "gas-impermeable" membranes have been known for 
almost 25 years, it is only 15 years ago that the first reports sug- 
gesting the existence of protein gas channels have appeared. In the 
first part of this chapter we will discuss studies that have been in 
favor of the function of AQP1 and other proteins as channels 
for C0 2 : 

(1) Nakhoul et al. (1998) reported that expression of aquaporinl 
(AQP1) in Xenopus oocytes increases the rate of intracellular 
acidification seen upon exposure of the oocyte to an increase 
in extracellular pco 2 ■ They proposed that AQP 1 acts as a CO2 
channel. 

(2) Forster et al. (1998) in the same year reported that the effect 
of DIDS (4,4'- diisothiocyano-2,2'-stilbenedisulfonic acid) 
on the C 18 O le O-exchange kinetics in human red cells as 
observed by mass spectrometry cannot be explained solely 
by inhibition of the anion exchanger 1 . When applying the 
mathematical model that describes this process in terms not 
only of CA activity and HCO^ permeability, but also of CO2 
permeability, as first introduced by Wunder and Gros (1997, 
1998) and Wunder et al. (1998), they concluded that DIDS is 
a powerful inhibitor of red cell membrane CO2 permeability. 
Since DIDS is an E-amino group reagent affecting a variety 
of proteins, they concluded that CO2 permeation in the red 
cell membrane must be a process mediated by a membrane 
protein. 

(3) Again in the same year, Prasad et al. (1998) reported that 
reconstitution of AQP1 into liposomes increased their CO2 
permeability as determined from stopped-flow experiments. 
However, as discussed above, the absolute values of the CO2 
permeabilities these authors obtained (lCP 3 cm/s) suggests 
that technical problems with the stopped-flow technique 
caused a problem in these measurements. Actually, the same 
group later on argued very clearly against a role of AQP1 as 
a channel for C0 2 (Missner et al, 2008; Zocher et al., 2012). 
So these results should actually not be used as an argument 
in favor of gas channels. 

(4) An important argument in favor of aquaporin's role as a CO2 
channel came more recently from the work of Uehlein et al. 
(2003). They showed that the tobacco aquaporin NtAQPl, 
like AQP1 (Nakhoul et al., 1998), increases CO2 permeation 
when expressed in Xenopus oocytes. Moreover, NtAQPl was 
shown in plants to increase both cellular water and CO2 per- 
meability and promote leaf growth. Thus, the authors showed 
in plants for the first time a clear physiological function of a 
protein C0 2 channel. 

In the years 2006-2008 Endeward et al. (2006a,b, 2008) elu- 
cidated the mechanism of CO2 permeation across the human 



red cell membrane. In these studies, using blood from AQP1- 
and Pvh-deficient persons, they were the first to detect a second 
C0 2 channel in addition to AQP1, namely the PJiesus pro- 
tein PJiAG (Rhesus-associated glycoprotein), which was already 
known to conduct ammonia (Ripoche et al., 2004). AQP1 and 
RhAG each contribute one half of the normal C0 2 permeabil- 
ity of this membrane, which amounts to Pco 2 = 0.15 cm/s, about 
10 times greater than the permeability in the functional absence 
of the two channels. Both channels are efficiently though not 
quite completely inhibited by DIDS, explaining the strong effect 
of this inhibitor on red cell Pco 2 an d thus the earlier finding by 
Forster et al. (1998). The C0 2 conductance through AQP1 was 
in addition strongly inhibitable by the mercury reagent pCMBS 
(p-chloromercuribenzene sulfonate) that is known to effectively 
inhibit the aquaporin water channel (Echevarria et al., 1993). In 
the Rh protein mediating CO2 transport, Endeward et al. (2008) 
demonstrated a competition between NH3 and CO2, which was a 
strong argument for RhAG being the C0 2 channel, rather than 
other members of the Rhesus protein complex in the red cell 
membrane. 

The evidence in favor of C0 2 channels listed so far has mostly 
been obtained on specialized membranes such as those of red cells 
and oocytes. Therefore, it was of interest to extend this to other 
membranes. Itel et al. (2012) have recently expanded the range of 
membranes exhibiting C0 2 conductance by aquaporin 1 by (a) 
incorporating AQP1 into artificial lipid bilayers containing 50% 
cholesterol and (b) expressing AQP1 in MDCK cells that contain 
about 40% cholesterol (see discussion below). 

The role of AQP1 as a CO2 channel was supported by molec- 
ular dynamics simulations by Hub and de Groot (2006) and by 
Wang et al. (2007), who both predicted that AQP1 lets pass molec- 
ular C0 2 , through the water channel of the AQP1 monomer, and 
perhaps mainly through the central pore that exists in the AQP1 
tetramer between the four monomers. The authors' reservation 
against a possible physiological significance of the CO2 channel 
came from their observation that the lipid bilayer surrounding 
the AQP tetramers was even more permeable to CO2 than AQP. 
Very recently, Uehlein et al. (2012) have provided an especially 
convincing piece of experimental evidence for AQP being indeed 
a channel for CO2. They incorporated the tobacco NtAQPl into a 
planar layer of gas-tight triblock-copolymers. While the untreated 
membrane shows no CO2 permeation whatsoever, incorporation 
of the NtAQPl induces a most impressive permeability for C0 2 . 
Hub et al. (2010) confirmed by molecular dynamics calculations 
that the Rhesus protein Rh50, a member of the Amt/MEP/Rh 
family, indeed can serve as a pathway for C0 2 in addition to NH3, 
thus confirming the experimental results on RhAG by Endeward 
etal. (2006b, 2008). 

Musa-Aziz et al. (2009) and Geyer et al. (2013a) showed, 
by expression in oocytes, that several members of the aqua- 
porin family 0-9 besides AQP1 are good conductors of CO2. The 
authors take the amplitudes of the (alkaline) surface pH tran- 
sients ApHs that appear upon exposure of the oocytes to CO2 as 
a semiquantitative measure of the aquaporin-induced CO2 per- 
meability. Somersalo et al. (2012) presented an extensive theory 
deriving the Pco 2 of the oocyte from the value of the transient 
ApHs peak. Somersalo et al. (2012) point out the limitations of 
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this measure of Pco 2 to be due to (a) the the peak value of the pH 
surface transient not being linearly related to Pco 2 > (b) the blunt 
surface pH electrode being pushed up tightly against the vitelline 
membrane of the oocyte, which creates a limited unstirred space 
there (of a size not precisely known) and causes the pH transient 
sensed by the electrode to be likely different from what it is in 
the remainder of the oocyte surface, and (c) the presence of the 
vitelline membrane per se being expected to create an additional 
unstirred space on the oocyte cell membrane, whose influence is 
hard to assess. Thus, comparisons of these pH transients caused 
by different aquaporins expressed in oocytes will be possible and 
provide most valuable insights, but will be semiquantative rather 
than quantitative. With this approach, the mentioned group finds 
ApHs to be largest for AQP5, and just slightly less for AQP6 and 
AQP1, such that these three seem to be the best conductors of 
CO2. Substantial ApHs values are also exhibited by AQPO, AQP4- 
M23 (a variant of AQP4-M1, which forms orthogonal arrays of 
tetramers in membranes), AQP9, and perhaps AQP8. Little or 
no C0 2 permeability is contributed by AQP2, AQP3, AQP4-M1, 
and AQP7. Functionally, it is notable that high permeabilities for 
CO2 seem to be exhibited by the most widespread isoform AQP1, 
occurring for example in red cells and the kidney, by the AQP4- 
M23, which is present in high density in the astrocytic endfeet 
at the blood-brain barrier, by AQP5, which is strongly expressed 
in alveolar type I pneumocytes and thus may facilitate diffusion 
of CO2 across the alveolar-capillary barrier, and finally by AQP6, 
which is present and colocalized with H + -ATPase in the intra- 
cellular vesicle membranes in the renal collecting duct epithelia. 
According to the authors' speculation, AQP6 may serve there to 
mediate CO2 exit from the vesicle lumen, after the latter has 
actively taken up H + , followed by its reaction with HCO^ to 
form CO2 inside the vesicle. Recently, Geyer et al. (2013b) have 
shown with their surface pH transient technique that in addi- 
tion to RhAG the tissue Rh proteins RhBG and RhCG serve as 
pathways for CO2. 

In the following second part of this chapter we will present 
some work that has come to the conclusion that either CO2 
channels do not exist or are functionally insignificant: 

Yang et al. (2000) studied the question of the function of 
AQP1 as a potential CO2 channel in a variety of systems, and 
reached a negative answer for all of them. They studied lipo- 
somes into which AQP1 was reconstituted, examined them by 
stopped-flow spectrophotometry, obtained the very low Pco 2 a s 
already discussed above and observed no effect of aquaporin- 
1 incorporation. Since the low Pco 2 suggests a major unstirred 
layer or mixing problem, as discussed, it appears possible that 
an action of AQP1 may have been obscured by the dominating 
effect of unstirred layers. The same holds for the authors' deter- 
mination of Pco 2 of red cells from normal and AQPl-ko mice, 
which was also done by stopped-flow. Their values of 0.012 and 
0.01 1 cm/s are more than 10 times lower than the values obtained 
by ls O mass spectrometry (Endeward et al., 2006a, 2008) and 
also than earlier rough estimates of red cell Pco 2 by Roughton 
(1959) and Forster (1969), yielding 0.15 and 0.58 cm/s, respec- 
tively. Again, it appears quite possible that Yang's low values are 
due to an unstirred layer or mixing problem in the stopped-flow 
apparatus, which would render changes in true membrane Pco 2 



insignificant in comparison to a predominant diffusion resistance 
in the extracellular solution. A third line of evidence presented 
by Yang et al. (2000) consists of an experiment with anesthetized 
artificially ventilated mice. These animals were subjected to a step 
decrease in inspired pco 2 from 35 to 0 mmHg. The absolute value 
and kinetics of the ensuing decrease in arterial pco 2 was then 
followed in normal and AQPl-ko mice, and no differences in 
either pco 2 amplitude or kinetics were observed. To appreciate the 
implications of this result, one has to recall the steps that could 
be crucial for the process of CO2 equilibration between blood 
and red cells. Firstly, red blood cells are fundamentally involved 
in this equilibration possess by chemical reactions and transport 
processes, and they possess AQP1 in their membrane. However, 
as discussed by Endeward et al. (2008), the mere lack of AQP1 in 
the red cell membrane (with RhAG still present) reduces Pco 2 to 
about 0.07 cm/s, and these authors showed that this value under 
resting conditions still allows complete CO2 equilibration of the 
red cell within the capillary transit time. Thus, red cells of AQP1 
null mice are not expected to behave functionally differently from 
normal mice in a gas exchange experiment under resting condi- 
tions. The second potentially crucial step is the transfer of CO2 
across the alveolar-capillary barrier. Among the elements of this 
barrier, only the endothelial membranes are candidates, as they 
express AQP1, while the alveolar epithelium does not (Verkman, 
1998; Matsuzaki et al., 2009). The latter, as mentioned above, 
exhibits a strong expression of AQP5. We will analyse the chances 
of observing an effect of the AQP1 deficiency in the alveolar- 
capillary barrier in the experiment of Yang et al. (2000) by using 
the classical concept of the lung diffusing capacity that goes back 
to Roughton and Forster (1957). Accordingly, we separate the 
resistance to CO2 transfer from lung to red cell into a "diffusion, 
or membrane" component and a "chemical reaction" component. 
This is expressed by I/Dl = I/Dl.m + 1/9 Vc, where Dl is the 
overall diffusing capacity of the lung, the reciprocal of Di t m is 
the diffusion or "membrane" resistance comprising that of the 
alveolar-capillary barrier, 0 is the chemical reaction rate in the 
red cell, and Vc the volume of the lung capillary blood. In other 
words, l/0Vc is the resistance to CO2 transfer provided by the 
chemical reactions of CO2 and associated erythrocytic transport 
processes. The comparatively small contribution of the diffusion 
resistance across the alveolar-capillary barrier becomes manifest 
in the minimum estimate for -Di m, HOOml/mmHg/min, which 
compares with the overall diffusing capacity of the lung Dl of 
180 ml/mmHg/min (Schuster, 1987). Thus, Di m makes up at 
most 1/6 of the total resistance for CO2 exchange between lung 
and red cells. In addition, the endothelial membranes make up 
only a fraction of the diffusion resistance of the alveolar-capillary 
barrier, 1 /Di t m- We conclude from this that, in order to become 
noticeable in measurements of CO2 exchange between blood and 
alveolar gas, loss of AQP1 from the endothelial membranes would 
have to cause a dramatic increase in the diffusion resistance of 
these membranes. Such an enormous effect of AQP1 seems not 
very likely. The same conclusion has been reached by Swenson 
et al. (2002), who attempted to inhibit pulmonary AQP by Zn ++ 
and, seeing no effect on CO2 exchange between diluted blood 
(10% hematocrit) or red cell-free lung perfusate and alveolar gas, 
estimated that the alveolar-capillary membrane diffusing capacity 
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Dl, m would have to be reduced by >90% to lower CO2 exchange 
by 10%. These arguments also apply to the observation of Fang 
et al. (2002), who in mice lacking AQP1 and/or AQP5 saw no 
difference to normal mice when they studied the exchange of 
CO2 between fluid-filled alveoli and lung-vessel perfusates. In 
that experiment, although of great interest because it studies lungs 
without AQP1 as well as AQP5, it appears likely that the large 
unstirred fluid zones in the alveoli will drastically decrease appar- 
ent m- Therefore, these unstirred fluid zones will probably 
dominate the process of CO2 exchange rather than the perme- 
ability of the membranes of the alveolar-capillary barrier. We 
conclude that the CO2 exchange experiments conducted in lungs 
do not seem to be quite conclusive evidence against the function 
of AQP1 (or AQP5) as CO2 channels, although the quantitative 
physiological role of these channels in gas exchange in the lung 
clearly remains debatable and should be studied further. 

A short note on the role of AQP1 in the kidney: Fang et al. 
(2002) have isolated proximal tubular vesicles, which are known 
to strongly express AQP1, and measured their CO2 permeability 
by the stopped-flow technique. As for the other vesicles studied by 
this group and already discussed above, the Pco 2 they obtained 
was very low, 0.0035 cm/s, and it was identical between vesicles 
prepared from normal and AQP1 -deficient mice. So the reser- 
vations against the stopped-flow technique explained above also 
apply to these experiments. On the other hand, Boron (2010) 
presents a preliminary report by Zhou et al. (2006), describing 
experimental work on intact proximal tubules. They come to the 
conclusion that 60% of the CO2 absorbed in the proximal tubule 
seems to be due to AQP1. 

Missner et al. (2008) in a more recent paper presented essen- 
tially two arguments against a functional role of CO2 channels 
in membranes. Firstly, they reported, as mentioned above, a CO2 
permeability of 3.2 cm/s in planar lipid bilayers, measured at an 
extremely alkaline pH of 9.6. At this alkaline pH, in the presence 
of CO2, HCO3" and CO3"-, CO2 diffusion within the considerable 
unstirred layer of >100|xm was accelerated due to a facilitated 
diffusion mechanism that had also been utilized for such mea- 
surements by Gutknecht et al. (1977). This is expected to result, 
as the authors illustrate theoretically, in an only modest pco 2 g ra ~ 
dient in the unstirred layer adjacent to the planar lipid bilayer, 
but in a significant pco 2 drop across the bilayer. This drop in 
Pco 2 is an expression of the transmembrane resistance for CO2. 
It also is the basis of the figure of 3.2 cm/s they obtain for Pqo 2 
from pH profile measurements in the front unstirred layer. The 
experimental procedure per se was similar to the one described 
above for the experiments of Zocher et al. (2012). We point out, 
however, the limitations of the fitting procedure used by Missner 
et al. (2008) to derive Pco 2 values, as already voiced in chapter I 
of this article. Missner et al. (2008) take the Pco 2 value of 3.2 cm/s 
to be representative for cell membranes in general, and go on to 
argue that with an unstirred layer (USL) of 1 \im as one often has 
around cells, the resistance of the USL vs. that of the membrane 
will be: 

d\jsh/D co 2 + l/Pco 2 = total apparent membrane diffusion 
resistance, or, inserting numbers: (7a) 



1 • 10 -4 cm/1.8 • 10~ 5 cm 2 /s+ 1/3.2 cm/s 

= 5.6s/cm + 0.3s/cm = 5.9s/cm. (7b) 

This shows, according to the authors, that the resistance to CO2 
diffusion through a cell membrane will always be dominated by 
the unstirred layer effect. Therefore, they conclude, a CO2 chan- 
nel in the cell membrane can never be of functional significance. 
We note that the basis of this argument rests on the assumption 
that the -Pco 2 °f the cell is identical to the Pco 2 of a lipid bilayer. 
This is in fundamental contrast to the determinations of red cell 
Pqo 2 (Endeward et al., 2006a, 2008), where even in the presence of 
CO2 channels values of 0. 1 5 cm/s were found, and in their absence 
an estimate of ~0.01 cm/s was obtained. It is also in contrast to 
the findings of Itel et al. (2012) who report a dramatic decrease 
of bilayer Pco 2 over 2-3 orders of magnitude with increasing 
cholesterol content, and it is in contrast to the Pco 2 of 0.01 cm/s 
observed for cells in culture (Itel et al., 2012) as well as to earlier 
measurements of a Pco 2 of 0.001 cm/s for the apical membrane of 
colon epithelium (Endeward and Gros, 2005). The calculation of 
Equation (7a), when performed for a Pco 2 of 0.01 cm/s, leads to a 
25 times greater diffusion resistance of the membrane in compar- 
ison to that of the 1 |im thick unstirred layer. Thus, the available 
experimental Pco 2 values of cell membranes suggest that protein 
CO2 channels can indeed be physiologically significant. 

The second argument presented by Missner et al. (2008) 
against the relevance of CO2 channels in biological membranes 
comes from measurements of the flux of CO2 across mono- 
layers of MDCK cells. They study native MDCK cell layers in 
comparison to MDCK cells expressing AQP1 and substantiate 
the presence of AQP1 by immunocytochemistry and by water 
flux measurements. In terms of CO2 flux, they find no differ- 
ence between MDCK cells with and without AQP1. In order to 
appreciate this latter result, one has to analyse the experimen- 
tal situation employed by Missner et al. (2008). We do this in 
the following along the lines already published by Endeward and 
Gros (2009). The cells are grown on semipermeable supports by 
Transwell (Corning), whose exact type is not given, but which 
may have been 50 u,m thick. The flux measurement then was 
done across the sum of the MDCK monolayer, which may have 
been 10 \im thick, plus the 50 u,m support, plus an unstirred layer 
on the back side, behind the support, of the thickness demon- 
strated by the authors for the front side to be at least 130 |xm. 
This adds up to a total water phase of at least 190 |xm in thickness, 
across which the flux of CO2 occurred. We ignore the front side 
unstirred layer, because that was controlled by the microelectrode 
measurement and gave the authors the precise front-side surface 
pH value. If we now perform, as above in Equation (7), an anal- 
ysis of the diffusion resistances of the overall water phase and of 
the membranes, we can write: 

d w /Dco 2 + l/Pco 2 b + l/Pco 2 a = total diffusion resistance 
of the combined layer, (8) 

where d w is the overall thickness of the combined layers, Dco 2 
the CO2 diffusion coefficient in water (taken to approximate dif- 
fusion in all three components of the layer), -Pco 2 b the CO2 
permeability of the basal MDCK membrane, and Pco 2 a the CO2 
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permeability of the apical membrane. To insert numbers, let us 
first use Missner's assumed Pco 2 f° r cells, 3.2 cm/s: 

190- 1(T 4 cm/1.5 - l(T 5 cm 2 /s+ 1/3.2 cm/s + 1/3.2 cm/s 
= 1270 s/cm + 0.3 s/cm + 0.3 s/cm = 1271 s/cm (8a) 

This would imply that the membrane resistance makes 
up <0.05% of the total resistance, and no change in membrane 
Pco 2 by a factor even of 10 would ever be detectable. Thus, if we 
accept the logic of the paper by Missner et al. (2008), an aqua- 
porin effect on membrane permeability could never have been 
measurable. Thus, the conclusion that AQP1 does not increase 
membrane CO2 permeability is not substantiated by this experi- 
ment. 

Instead of Pco 2 = 3-2 cm/s, we can tentatively insert into 
Equation (8) the value determined by Itel et al. (2012) for isolated 
MDCK cells, P C o 2 = 0.017 cm/s: 

1270 s/cm + 1/0.017 s/cm + l/0.017s/cm = 1390s/cm(8b) 

Even with this much lower membrane permeability, the two 
membranes contribute < 10% of the total diffusion resistance. An 
increase in Pco 2 due to AQP1 expression by ~50% as reported 
by Itel et al. (2012), would necessitate being able to measure the 
difference between total diffusion resistances of 1390 s/cm and 
1350 s/cm, i.e., a difference of less than 3%. In view of the preci- 
sion of these measurements — and of any other conceivable exper- 
imental technique — it seems impossible to detect such an effect of 
AQP1 on membrane Pco 2 ■ No matter which membrane Pco 2 we 
accept, the experiment used by Missner et al. (2008) seems to be 
unsuitable to observe AQP effects on Pco 2 ■ The decisive flaw in 
this setup is the large overall thickness of the diffusion layer. 

Summarizing this chapter with the pros and contras on the 
existence of CO2 channels, we feel that overwhelming evidence 



has been accumulated by now that aquaporins and PJiesus pro- 
tein indeed act as CO2 channels, and several lines of evidence 
proposed against CO2 channels are subject to some critique. 
Whether these channels are physiologically important or not, 
indeed depends crucially on the question what the "intrinsic" 
CO2 permeability of a cell membrane is. If it is very high, chan- 
nels do not make sense, if it is low, they do. A small paper 
by representatives of the opposing views on this topic, Boron, 
Endeward, Gros, and Pohl, based on their discussion at a meet- 
ing in Strobl (Austria), was published (Boron et al., 2011) and 
summarizes this state of the discussion. After this meeting, the 
paper by Itel et al. (2012) has been published that shows that 
many membranes, due to their cholesterol content, have rather 
low intrinsic CO2 permeabilities and thus channels do make 
sense. This paper demonstrates this by (a) incorporating AQP1 
into artificial phospholipid vesicles containing 50% cholesterol, 
which increases the vesicles' Pco 2 up to 9-fold (Figure 7A), and 
(b) by showing that expression of AQP1 in MDCK cells, which 
possess ~40% cholesterol, significantly increases membrane CO2 
permeability by ~50% (Figure 7B). Both experiments clearly 
illustrate that AQP1 acts as a functionally relevant CO2 channel, if 
the "intrinsic" membrane permeabilities for CO2 are sufficiently 
low, ~0.003 cm/s in the case of vesicles and 0.017 cm/s in the case 
of the MDCK cells. As most cells have cholesterol contents of 
around 40%, most cells are expected to exhibit an intrinsic CO2 
permeability in the order of 0.01 cm/s, and thus their CO2 per- 
meability will be significantly enhanced by the incorporation of a 
protein gas channel. It appears therefore at this stage that nature 
has two options to achieve a high CO2 permeability of a cell mem- 
brane, either by reducing its content of cholesterol, which will be 
at the expense of a reduced barrier function to other molecules 
and of a reduced mechanical stability, or by maintaining a high 
cholesterol content and raising CO2 permeability by incorporat- 
ing protein gas channels. It will be of interest in the future to find 
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FIGURE 7 I Effects of AQP1 incorporation into lipid vesicles and into 
MDCK cell membranes. (A) Phospholipid vesicles of 
phosphatidylcholine:phosphatidylserine 8:2, with 50 mol% cholesterol, 
with and without incorporation of human AQP1. LPR is lipid-to-protein 
ratio. At the highest protein incorporation, LPR 140, the P C o z as 
measured by 18 0 exchange mass spectrometry is 9-fold increased over 
the control (leftmost column). DIDS reduces the AQP1-mediated CO2 



conductances, but not the Pco 2 of control vesicles. ### Signif icantly 
different from AQP1-free control vesicles, P<0.01; 'significant effect of 
DIDS on vesicles with the same LPR, P < 0.05. (B) Expression of 
human AQP1 in MDCK cells increases P C o z significantly (**P<0.02). 
DIDS reduces the hAQP1-mediated increase in P C02 ( $S P<0.02), but 
does not affect Pco 2 of control cells. Reproduced from Itel et al. 
(2012), with permission. 
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more biological examples of these two options. A third option is 
of course to combine high cholesterol content with the absence of 
gas channels, which will give "gas-impermeable" membranes as 
discussed above. 

A BRIEF VIEW ON OTHER GASES OF PHYSIOLOGICAL 
INTEREST 

Direct determinations of cellular gas permeabilities for the gases 
listed in Table 1 other than CO2 are sparse. Kutchai (1975) con- 
cluded that the permeabilities of red cell membranes for O2 and 
CO are high enough in order to constitute no significant dif- 
fusion resistance during O2 and CO uptake by these cells, but 
did not come up with numbers for these permeabilities. Honig 
et al. (1992) observed by oxymyoglobin microspectrophotometry 
a shallow po 2 gradient inside skeletal muscle cells and concluded 
that a steep gradient of po 2 must occur between the surface of 
the red cell and the sarcolemma. In this phenomenon, however, 
the authors implied neither the red cell membrane nor the sar- 
colemma as sites of major diffusion resistance to O2, although 
this would have been conceivable, but rather attributed this resis- 
tance to the interstitial myoglobin-free space. However, in a later 
paper Voter and Gayeski (1995) found that the spatial resolu- 
tion in the earlier experiments was considerably poorer than 
previously thought, so the intracellular po 2 gradients were in 
fact greater, reducing the role of the space between red cell 
and muscle cell membrane as a diffusion resistance for O2. In 
the present context, these data do not give an indication for a 
significant role of the two membranes as diffusion resistances 
toward O2. For lipid vesicles, Widomska et al. (2007) used the 
saturation-recovery EPR technique to determine the oxygen per- 
meability coefficient. They reported extremely high Pq 2 values 
of 51, 50, and 157cm/s for lipid bilayers at 37° C made from 
lipids extracted from membranes of the lens, from a 1 -palmitoyl- 
2-oleoylphosphatidylcholine/cholesterol mixture, and from pure 
l-palmitoyl-2-oleoylphosphatidylcholine, respectively. The same 
group has reported measurements with this technique in the 
Chinese hamster ovary plasma membrane, finding a Po 2 °f 
42cm/s at 37°C (Subczynski et al, 1992). These values, if true 
for cell membranes in general, would be far beyond constituting 
a barrier to O2 supply to cells, and they would deviate drasti- 
cally from the value predicted in Table 1 from the measured Pco 2 
and the gas partition coefficients. If the very high intramem- 
brane diffusivity of O2 arising from these high Po 2 values were 
true for O2 and perhaps other gases such as NO, CO, and N2 as 
well, the implication would be that even in a pure phospholipid 
membrane, Dm for CO2 would be exceptionally low among these 
gases. In sharp contrast to these EPR measurements, Ivanov et al. 
(2004) reported for surface lipid monolayers oxygen permeabil- 
ities between 10~ 4 and 10^ 5 cm/s, i.e., 6-7 orders of magnitude 
lower than the values derived from EPR. As a consequence, they 
postulated the existence of O2 channels in biological membranes. 
Indeed, Wang et al. (2007) have shown by molecular dynamics 
simulations that AQP1 can serve as a conduit for O2. The group of 
Echevarria et al. (2007) has also postulated the existence of oxygen 
channels, not from direct O2 flux measurements, but based on 
several lines of indirect evidence, e.g., they report (a) that overex- 
pression of AQP 1 results in an accelerated upregulation of HIF2a 



after exposure to hypoxia of PC12 cells, (b) that knockdown of 
AQP1 in endothelial cells induces hypoxia-inducible genes, and 
finally (c) they observe that rat lungs under systemic hypoxia up- 
regulate the expression of AQP 1. One might ask, though, whether 
some of these effects could not be related to the function of AQP1 
as a water channel. 

A somewhat similar situation exists with respect to NO. Sakai 
et al. (2008) studied the uptake of NO by hemoglobin-loaded arti- 
ficial vesicles using stopped-flow spectrophotometry. With low 
intravesicular hemoglobin concentrations and small vesicle sizes, 
the rate of combination of NO with deoxyhemoglobin was iden- 
tical to the reaction rate with deoxyhemoglobin in solution. The 
authors concluded that no resistance of the vesicle membrane 
was apparent. With increasing vesicle size a retardation of the 
apparent association rate constant became apparent, which they 
interpreted as being due to an increasing role of intravesicular gas 
diffusion. This latter effect was not visible when they measured 
the uptake of CO by these vesicles, a fact they attributed to the 
much slower reaction kinetics of the CO-hemoglobin association. 
Again, in the case of NO, molecular dynamics simulations pre- 
dict that AQP4 is an especially good conductor of both NO and 
O2 (Wang and Tajkhorshid, 2010), which would not make sense 
physiologically if the membrane permeabilities were generally 
immeasurably high. Along this line, in the case of NO a substan- 
tial body of evidence has been compiled from AQPl-transfected 
cells in culture, from AQP1 reconstituted into lipid vesicles, and 
from physiological measurements of the AQP 1 -dependence of the 
physiological actions of NO, all indicating that AQP1 is essential 
in getting NO across cell membranes from its site of production 
to its target (Herrera et al, 2006; Herrera and Garvin, 2007). 

In conclusion, the field of O2 and NO permeation across 
biological membranes is highly controversial and drastically con- 
tradictory results have been reported. This area clearly awaits 
further studies of the mechanisms by which these gases cross 
membranes. 

A DISCUSSION OF THE SUITABILITY OF THE 18 0 EXCHANGE 
MASS SPECTROMETRIC TECHNIQUE TO MEASURE THE 
FAST PROCESS OF C0 2 PERMEATION 

We have discussed above that CO2 exchange with lipid vesicles of 
a diameter of 150 nm may not be measurable by rapid reaction 
stopped-flow techniques, which have dead-times of a few mil- 
liseconds, provided the mixing efficiency of the apparatus is very 
good. This is so because the process is much faster than the dead 
time, if intravesicular CO2 hydration velocity is not limiting and 
if there is no marked diffusion resistance offered by the vesicle 
membrane. How can under these premises the mass spectromet- 
ric ls O technique overcome the problem of the extreme rapidity 
of the processes of CO2 permeation? In this chapter we will first 
summarize the principle of the ls O exchange technique, which 
has already been described before in great detail, then address 
the above question how a mass spectrometric technique, whose 
response time is 3 s (Endeward and Gros, 2005) and which records 
its signal on a time scale of several hundreds of seconds, can be 
able to describe the extremely fast process of CO2 exchange across 
cells and vesicles. Next, we will demonstrate the sensitivity of the 
calculated permeability values for CO2 to other parameters that 
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FIGURE 8 | Schematic representation of the processes occurring when 
carbonic anhydrase-containing cells are added to 1s O-labeled 
CO2-HCOJ solution. The labeled O is marked in blue, the most critical 
parameters determining the time course of the process are labeled in red. 
Single arrows indicate transport processes, double arrows reaction 
processes. 
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FIGURE 9 I Original mass spectrometric recording of an experiment 
with human red cells. The concentration of C 18 0 16 0 is plotted vs. time. 
First arrow point of addition of red cells, second arrow point of addition of 
CA in excess to rapidly establish final isotopic equilibrium. Solid curve (red): 
mass spectrometric recording; dotted curve (yellow): calculated curve using 



the best-fit values of Pro, and P u , 



. 37C. 



enter into the calculation, then illustrate quality and reliability 
of the fitting procedure that needs to be employed, because an 
analytical solution of the system of differential equations describ- 
ing the ls O exchange process is not available, and finally discuss 
experimental data on the effect of unstirred layers in the mass 
spectrometric Pco 2 measurements of red cells. 

PRINCIPLE OF THE METHOD 

We will briefly summarize the mass spectrometric method to 
measure Pco 2 of cells and vesicles, which has first been described 
in detail by Wunder and Gros (1997, 1998) and Wunder et al. 
(1998), and later in its final form by Endeward and Gros (2005). 
All technical details are described in these papers. The underly- 
ing principle is that used by Mills and Urey (1940) to measure 
the kinetics of the CO2 hydration-dehydration reaction by using 
ls O-labeled CO2 and HCO^ and observing the very slow process 
of the disappearance of the ls O from the CO2/HCO3" pool into 
the comparatively huge H2O pool. Much later, Itada and Forster 
(1977) developed an ingenious gas inlet system for the mass 
spectrometer and used this method to observe the effect of the 
addition of red cells to the solution of labeled CGVHCO^upon 
the time course of ls O loss from CO2. They calculated from this 
the intraerythrocytic carbonic anhydrase activity and the red cells' 
bicarbonate permeability. Wunder and Gros (1997, 1998) and 
Wunder et al. (1998) became aware of the fact that this time 
course also contains information on the CO2 permeability of 
the red cell membrane and they developed the theory describ- 
ing this process, which allowed them to derive Pco 2 values from 
the mass spectrometric recordings. Because some refinements of 
the theory were added later on, we will in the following refer to 
the theory in the form described by Endeward and Gros (2005). 
NaHC 18 0 16 02 is dissolved in a solution in the reaction chamber 
connected to the mass spectrometer via the previously described 
gas inlet system (Wunder and Gros, 1998). The ls O-labeled bicar- 
bonate rapidly establishes chemical equilibrium with ls O-labeled 
CO2, and from then on chemical equilibrium exists in the sys- 
tem, and the entire remainder of the experiment occurs under 
conditions of perfect chemical equilibrium, including constant 
values of pH and total pco 2 ■ However, the system at this point 
is far from isotopic equilibrium, which is approached very slowly 
by the mechanism indicated in the upper half of Figure 8. Each 
time the labeled bicarbonate reacts to give water and CO2, there 
is a 1/3 chance that the ls O is lost from the CO2-HCO3 pool 
into the water pool. Since the latter is 55.5 M/20 • 10~ 3 M » 3000 
times greater than the former, eventually almost all the 18 O will 
be lost into the water. This causes a slow spontaneous decay of 
C 18 O le O, the quantity measured by the mass spectrometer via 
the gas inlet system. This decay is seen in the first part of the 
recording of Figure 9. The slope of this first part is what Mills 
and Urey (1940) used to determine the CO2 hydration velocity 
constant. Both, Figures 8, 9 show that the situation changes dra- 
matically, when red cells (or other carbonic anhydrase-containing 
cells or vesicles) are added. This is followed by a (rapid) influx of 
labeled CO2 into the red cells, where it is fed into the same reac- 
tion processes as shown for the extracellular space in the upper 
half of Figure 8. More slowly than C 18 O le O, labeled bicarbonate 
then also enters the red cells and is fed into these reactions. 



These two processes can be assigned — although not exclusively — 
to the two phases seen in the mass spectrometric recording in 
Figure 9 after red cell addition: the first rapid phase is domi- 
nated by the CO2 influx (besides by the intracellular carbonic 
anhydrase (CA) activity), and the second slower phase contains 
most of the information about the influx of HCO^ . The same 
reaction processes that occur in the extracellular space occur in 
the intracellular space 20,000 times faster due to the CA, and 
for this reason almost the entire reaction process now takes place 
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intracellularly. As evident from the scheme of Figure 8, the intra- 
cellular reaction process depends on CO2 and HCO^ moving into 
the cells, a process which is controlled by the permeabilities of the 
membrane for CO2, Pco 2 > an d for HCO^, f^co -- All me reac- 
tion and diffusion processes of Figure 8 are described by a set of 
six differential equations (Endeward and Gros, 2005), and their 
numerical solution is used to determine by a fitting procedure 
the best-fit values of these two parameters. It should be noted 
that the intracellular CA activity, which is also very important for 
the kinetics of the ls O exchange process, is determined indepen- 
dently from cell lysates under intracellular conditions and is not 
obtained from the fitting procedure. Determining CA activity in 
this way requires ensuring that two important parameters influ- 
encing specific CA activity equal those present inside cells, i.e., 
pH and [Cl _ ], and it presupposes that the CA activity measured 
in diluted cell lysate can safely be extrapolated to the intracellu- 
lar CA concentration. This latter assumption has previously been 
validated by the work of Donaldson and Quinn (1974). The per- 
fect superimposition of the original mass spectrometric recording 
in Figure 9 (solid red curve) and the theoretical curve (dotted 
yellow) calculated with the best- fit values of Pco 2 an d ^hco - * s 
achieved for almost all experiments. It illustrates that the theory 
provides an excellent description of the process of Figure 8. 

HOW CAN THE 18 0 EXCHANGE MASS SPECTROMETRIC METHOD 
MEASURE PROCESSES THAT OCCUR WITHIN SOME MILLISECONDS? 

The difference between the kinetics of cellular CO2 exchanges 
and the kinetics of the decay of the ls O-labeled CO2 reflecting 
these exchanges is huge. For example, let us consider the CO2 
uptake process by human red cells. This process is so fast that 
it is tricky to measure it by the stopped-flow technique, for the 
reasons already discussed, dead time of a few milliseconds, possi- 
ble imperfect mixing, and the special problem encountered when 
mixing cell or vesicle suspensions in this apparatus. For this rea- 
son, stopped-flow is not an ideal technique to measure the uptake 
of gases by cells, and we will not refer to such measurements 
here. Holland and Forster (1975) measured the uptake of CO2 by 
human red cells with a rapid reaction continuous flow apparatus, 
which does not have this problem, because the solutions with 



the cells continue to flow during observation. The authors deter- 
mined a half-time of this process of 13 ms, which we consider 
fairly reliable for the reasons just given. Let us compare this num- 
ber with the kinetics of the uptake of C 18 O le O by human red cells, 
as it is reflected in the extracellular decay of this species and illus- 
trated in Figure 9. Using the point of red cell addition, and the 
final isotopic equilibrium reached after the addition of excess free 
CA at the end of the experiment, we can estimate the half-time 
of the overall process to be ~140s, which is more than 10,000 
times the half-time of net uptake of CO2. What is the mecha- 
nism of this enormously slow representation of a very fast process 
in the C 18 0 16 0 decay signal? There are two major mechanisms 
contributing to this, which we shall discuss in the following. 

(1) Mills and Urey (1940), who developed the ls O exchange 
technique to measure the kinetics of the uncatalyzed CO2 hydra- 
tion reaction, explained that the disappearance of ls O from the 
labeled CO2/HCO3" into the water pool will, firstly, be only 1/3 
as fast as the net CO2 hydration reaction as it is illustrated in 
Figure 10A, because under chemical, but not isotopic, equilib- 
rium only 1/3 of the reaction steps producing H2O and CO2 from 
HCO3" will lead to as loss of ls O into the water pool (see upper 
half of Figure 8). Secondly, as most of the ls O label is located in 
HCO3" rather than in CO2, but each forward and backward reac- 
tion step in Figure 8 has to pass through CO2, the process will be 
slowed by the ratio [HCO^j/ICC^], which for the calculation of 
Figure 10B we have assumed to be 17, corresponding to a pH of 
7.3 at 37°C. 

In Figure 10B we have calculated the half time necessary to 
reach isotopic equilibrium after a step increase in C 18 O le O in 
solution, in the manner in which this is also done for the extra- 
cellular space in the equations describing the process of Figure 8. 
This yields a ti/2 of 250 s. In Figure 10A we have calculated the 
half-time of the net reaction CO2+ H 2 0 -o- HCO^~+ H+ at 37°C 
after a step increase in |>co 2 > giving a ti/2 of 5 s. In other words, 
the decay of C 18 0 16 0 occurs 50 times more slowly than the net 
CO2 hydration. The factor 50 results from the product 3 x 17, as 
predicted by Mills and Urey (1940). 

Having accounted for a factor of 50 out of a total slow-down 
factor of 10,000, we are left with a factor of 10,000/50 = 200. 
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FIGURE 10 I Comparison of the time courses of the net CO2 
hydration reaction after a step increase in [C0 2 ] (A), and the 
spontaneous decay of C 18 0 16 0 after a step increase in labeled C0 2 



and HCOj being in chemical equilibrium but not in isotopic 
equilibrium with water (B). Calculated for 37°C. Inserts indicate the 
half-times of the two reactions. 
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FIGURE 11 I Calculated mass spectrometric recordings for three 
different conditions of hematocrit and intracellular activity. Blue curve: 
human red cell hematocrit of 0.02% with intracellular activity Ai of 20,000. 
Green curve: hematocrit of 0.2% with Ai = 20,000. Red curve 
(superimposed to the blue curve): hematocrit of 0.2% with Ai = 2000 and 
10 times smaller surface-to-volume ratio S/V, similar to what is exhibited by 
MDCK cells etc. Y-axis: actual concentration of C 18 0 16 0 at a given time 
minus its concentration [C 18 0 16 0]co after isotopic equilibrium has been 
reached. 



8 
o 



0.1 







PhC03" 


= 0.0015 cm/s 




\ x 


^\ ~~~~~~ — ~ 














^002 


= 0.15 cm/s /^^'N^/ *" 








P c02 = 0.01 5 cm/s / \\ 








P c02 = 0.0015 cm/s 










A e = 4\ 


) 


100 


1 1 

200 300 


i 

400 



t(s) 



FIGURE 12 | Calculated time courses of mass spectrometric recordings 
from human red cells at a hematocrit of 0.02%. The red curves are 

calculated for a P COz of 0.15 cm/s ( — ), 0.01 5 cm/s ( ) and 0.001 5 cm/s 

( ). The green line ( ) indicates an extracellular CA activity Ae of 4, 

if free CA is added to the solution in the absence of cells. If the red cells 
studied in the three red curves were lysed, they would also cause an Ae of 
4. Denotation of axes identical to Figure 11. 



(2) An important point to consider here is that in the mass 
spectrometric measuring chamber one does not use packed red 
cells, but highly diluted red cell suspensions. By observing the 
decay of the extracellular C 18 0 16 0 concentration, one does not 
observe what is going on inside the red cell, but its reflection in the 
outside solution. Obviously, this signal will vary drastically with 
the hematocrit in the measuring chamber. However, the optimal 
hematocrit is not just 1/200 = 0.5%, but even >10 times less 
(Endeward et al., 2008). The extracellular solution in view of this 
very low hematocrit contains a huge reservoir of ls O, mainly in 
the form of HC 18 0 16 C>2 • Due to the high intracellular CA activ- 
ity, the reactions of Figure 8, however, occur almost exclusively 
inside the cell. This implies that all the label of the extracellular 
solution has to pass through the very small volume of the red cells 
on order to end up eventually in the water of the reaction cham- 
ber. This situation will slow down greatly the overall process, and 
is comparable to an RC element with a capacitor of large capaci- 
tance (the large pool of ls O in the CO2/HCO3" of the extracellular 
solution) in conjunction with a resistor (the small red cell volume 
through which the process of ls O decay has to pass) in an electric 
circuit. Given that the time constant of an RC element is R ■ C, 
an increase in R will slow down the process. Inversely, decreasing 
R by increasing the hematocrit will accelerate the process. This 
latter case is illustrated in the curves of Figure 11, which were cal- 
culated with our mathematical model describing the scheme of 
Figure 8. The blue curve illustrates the often used condition of 
a red cell experiment with a hematocrit in the measuring cham- 
ber of 0.02% (with the usual intracellular CA activity of human 
red cells of 20,000, and a surface-to-volume ratio of these cells 
of 20,000 cm -1 ). The faster, green curve, which would just be 
measurable in the mass spectrometer, is calculated for a ten-fold 
hematocrit of 0.2%. It is apparent that this causes an enormous 
acceleration of the C 18 O le O kinetics, although not quite by a fac- 
tor of ten, but by a factor of about 5.5. This non-proportional 
behavior is the reason why, in order to achieve a deceleration of a 
factor of 200, one needs a greater than 200-fold dilution of packed 
red cells. 

In conclusion, the present mass spectrometric technique is 
able to measure extremely fast CO2 exchanges (i) because it 
observes the loss of 18 0 from CO2+ HCO^ into the water pool, 
rather than observing directly the CO2 exchange process, and (ii) 
because an extreme dilution of the cell suspensions drastically 
slows down this kinetics. An additional and practically highly 
useful information can be taken from Figure 1 1 by noting that 
the red curve coincides with the blue curve: cells with proper- 
ties approaching those of MDCK cells that are larger than red 
cells (with a correspondingly smaller surface-to-volume ratio), 
but possess, for example, a 10 times lower intracellular CA activ- 
ity, can be measured equally well, if in compensation the cytocrit 
is raised ten-fold. This implies that the method is suitable to mea- 
sure P C o 2 of cells like MDCK, HEK293 or tsA201 cells with con- 
siderably lower intracellular CA activity than red cells (Al-Samir 
et al., 2013), provided the cytocrit in the measuring chamber is 
raised correspondingly. 

The effect of using the extremely low hematocrits in these 
experiments is further visualized in Figure 12. With Ai = 20,000, 
a hematocrit of 0.02% implies an overall CA activity in the 
measuring chamber of 20,000 x 0.0002 = 4. The green line in 



Figure 12 shows the effect of adding an amount of free CA in 
solution into the chamber that yields a four-fold acceleration of 
the slope of C 18 0 16 0 decay. The red curves have been calculated 
for a hematocrit of intact red cells of 0.02% and an Ai of 20,000, 
using different Pco 2 values of 0.15, 0.015, and 0.0015 cm/s. The 
three curves are distinctly different, especially in terms of the 
kinetics and amplitude of their first phases, and thus conspicu- 
ously reflect the different Pqo 2 values. The overall CA activity of 
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the red cells studied in the three red curves, however, when aver- 
aged over the entire chamber volume, is 4. This activity when 
present free in solution is indicated by the green line (Ae = 4). 
It may be seen that the gross overall time course of all three red 
cell curves is similar to that of the green curve for Ae = 4. We 
conclude that the extreme dilution of the cells investigated, which 
this method allows us employ, makes it possible to slow the pro- 
cess of ls O exchange sufficiently down to make it measurable by 
mass spectrometry. 

SENSITIVITY OF THE P C o 2 VALUES ESTIMATED FROM THE MASS 
SPECTROMETRIC RECORDS FOR THE PARAMETERS USED IN THE 
THEORETICAL TREATMENT 

In addition to the two parameters to be determined, Pco 2 an d 
P HCO -, the parameters given on the x-axis of Figure 13 are used 
in the set of differential equations (Endeward and Gros, 2005). It 
is apparent that two very sensitive parameters are the intracellu- 
lar CA activity, Ai, and the extracellular pH, pHe, as it prevails 
in the mass spectrometric measuring chamber. As mentioned, Ai 
is carefully determined independently. pHe is equally carefully 
controlled by a pH electrode present in the measuring cham- 
ber and always kept precisely at pH 7.40 (Endeward and Gros, 
2005; Endeward et al., 2006a). K' v the first apparent dissociation 
constant of carbonic acid, is a precisely known physico-chemical 
constant. pHi, the intracellular pH, is usually well defined by the 
pHe and is entirely uncritical, as is the water permeability of the 
membrane, Ph 2 o- Even variations of Ph 2 o with AQP1 expression 
or lack thereof do not influence the results, because Ph 2 o is always 
high enough to warrant a sufficient efflux of labeled water from 
the cells, a process illustrated Figure 8. An important parame- 
ter next to Ai and pHe is the cellular surface-to-volume ratio a, 



which is determined on each experimental day for all cell types 
employed. In the case of red cells the volume is derived from 
hematocrit and cell density, and literature values are used for the 
surface area. In the case of isolated primary cells or cells isolated 
from cultures, we determine microscopically cell density and cell 
diameter and utilize their usually spherical shape to calculate a. 
In the latter case, the uncertainty is greater than in the case of red 
cells, but is not expected to exceed 10%, a number which would 
introduce an error of equal percentage into Pco 2 ■ 

In the context of discussing the influence of Ai on the determi- 
nation of Pco 2 > it seems appropriate here to add a note on how the 
ls O method behaves when cells of larger size but with much lower 
intracellular CA activity than red cells, such as MDCK cells, are to 
be measured. We consider cells with an intracellular CA activity 
of 500, i.e., 1/40 of that of red cells, which are employed at a suit- 
able cytocrit and possess a true Pco 2 °f 0-01 cm/s. We compare 
two curves of the kind shown in Figure 9, one fitted perfectly by 
Pco 2 = 0.01 cm/s, the other calculated with identical parameters 
except a value of Pco 2 of 0.1 cm/s forced upon the calculation. 
We determine from these curves the concentrations [C 18 O le O] 
shortly after the beginning of the 2nd slower phase, which fol- 
lows the 1st phase initiated by the addition of cells, and find 
that these concentrations differ between the two curves by 20%. 
This is a highly significant difference, which shows that in these 
cells, similar to the case of red cells (see Figure 1 in Endeward 
et al., 2008), CO2 permeabilities between 0.01 and 0.1 cm/s can 
be distinguished with great certainty. If we do the analogous cal- 
culations for Pco 2 =0.1 cm/s and Pco 2 = 1 cm/s, the difference 
in [C 18 0 16 0] is 3%, i.e., a distinction between these permeabili- 
ties is still possible, but with decreasing reliability. This behavior 
is similar to that exhibited by red cells. Testing the sensitivity of 
Pco 2 (0.01 cm/s) of these MDCK-type cells for changes in Ai, 
in the manner shown in Figure 13 for red cells, we find that an 
increase in Ai from 500 to 525 (+5%) causes a decrease in Pco 2 
by 11%, and a decrease in Ai from 500 to 475 (—5%) causes an 
increase in Pco 2 by 26%. Thus, for cells of low CA activity the 
sensitivity of Pco 2 for errors in Ai is equal to or lower than in 
the case of red cells as shown in Figure 13. The overall conclusion 
from these considerations is that the mass spectrometric method 
works equally well for cells with rather low Ai, e.g., of the type of 
MDCK cells, as it does for red cells. 

UNIQUENESS AND PRECISION OF THE PERMEABILITIES OBTAINED 
FROM THE FITTING PROCEDURE 

It is apparent from Figure 14A that the fitting procedure yields 
highly reliably the values of Pco 2 as well as P HCO - from a given 
experimental record. This is due to the above-mentioned fact that 
the curve of C 18 0 16 0 decay that evolves after the addition of the 
cells, exhibits two clearly separated phases. Each of the two per- 
meabilities does not solely determine one of these phases, but 
Pco 2 dominates the first rapid phase (together with Ai), while 
P HCO - plays a minor role, and P HCO - is dominant in the second 
slower phase, where Ai and Pco 2 play minor roles. This situation 
is extremely useful in unambiguously deriving the two permeabil- 
ities from the mass spectrometric records, and for the same reason 
the fitting procedure converges exceptionally well. Figure 14A 
shows that there is a unique minimum in the sum of squared 
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FIGURE 13 I Sensitivity of the estimated Pco 2 for variations in the 
parameters used in the system of equations decribing the processes 
of Figure 8. a, Surface-to-volume ratio of cells; K', , first apparent 
disscociation constant of carbonic acid: K!, = [H+] [HCO3 l/lCCy; Ai, 
intracellular CA activity; pHe, extracellular pH in measuring chamber; pHi, 
intracellular pH; Ph 2 o. water permeability of cell membrane. Reproduced 
from Endeward et al. (2008), with permission. 
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FIGURE 14 | Quality of the fitting procedure, expressed as sum of 
squared deviations between experimental and calculated 
concentrations of C 18 0 16 0. (A) Two-dimensional representation 
illustrating that a unique minimum of the sum of squares defines Pco 2 
as well as P H ccr f° r a given experimental record of C 18 0 16 0 decay. 
There are no local minima whatsoever. Reproduced from Al-Samir et al. 



(2013), with permission. (B) Sum of squares for an experimental curve 
with normal human red blood cells (blue curve) and with AQP1-deficient 
human red cells (red curve) as a function of Pco z - It is apparent that 
the different CO2 permeabilities of the two red cell types are well 
defined. The precision achieved by the fitting procedure for P C02 is 
better than 0.005 cm/s. 




Viscosity (10^ nf/s) 



FIGURE 15 | Plot of the reciprocal of the experimentally determined 



red cell P, 



co 2 



Cco 2 



) vs. the kinematic viscosity of the dextran 



solution in which it was measured. Extrapolation of the linear relation to 
zero viscosity yields the true membrane P C o 2 (Pco 2 tru0 )' 37°C. "Saline" 
indicates the Pco 2 measured in the absence of dextran. Reproduced from 

Endeward and Gros (2009), with permission. 



deviations between experimental and calculated concentrations 
of labeled CO2 for both Pqo 2 an d ^hco - > an ^ no l° ca l min- 
ima are present. Figure 14B illustrates for the examples of normal 
human red cells (blue curve) and of AQP1 -deficient human red 
cells (red curve) how well this method differentiates between the 
two cell types. Normal red cells in this example have a Pco 2 of 
about 0.1 1 cm/s, and AQP1 -deficient red cells of 0.07 cm/s, both 
values being in the range of values reported by Endeward et al. 
(2006a). The precision of these determinations is visibly better 
than 0.005 cm/s. 

ROLE OF UNSTIRRED LAYERS IN MASS SPECTROMETRIC 
MEASUREMENTS OF RED CELL C0 2 PERMEABILITY 

Because the diffusion resistance of cell membranes to CO2 
is in many cases — as discussed above — small, e.g., compared 
to that of ions, an additional resistance due to an unstirred 
layer of water around the cells studied is of importance. It 
is has therefore often been postulated that CO2 permeabili- 
ties were underestimated due to the presence of extracellular 
unstirred layers. Endeward and Gros (2009) have therefore devel- 
oped an approach allowing them to determine the unstirred 
layer thickness around red cells in mass spectrometric Pco 2 
determinations. 

From the additivity of the diffusion resistances of the mem- 
brane and the extracellular unstirred layer it follows that 

l/-Pco 2 , app = lAPco 2 , tme + S/^co 2 , (9) 

where Pqo 2 app i s the measured apparent CO2 permeability that 
includes effects by an unstirred layer, Pco 2 am is the true CO2 per- 
meability of the membrane, 8 is the thickness of the extracellular 
unstirred layer, and Dco 2 is the CO2 diffusion coefficient in water 
that describes CO2 diffusion in the unstirred layer. 

It is known from theoretical hydrodynamics that at given 
other conditions 8 increases with increasing solution viscos- 
ity v (Landau and Lifschitz, 1991). To utilize this relationship, 
Endeward and Gros (2009) suspended human red cells in solu- 
tions with dextran concentrations (mol.wt. 60,000) between 0 



and 10 g%, whose viscosities varied between 0.7 and about 4 • 
10~ 6 m 2 /s. At various dextran concentrations within this ranee 
they determined Pqo 2 app at 37°C, and plotted Pco 2 app vs - the vis- 
cosity of the dextran solution. The results are shown in Figure 15, 
which illustrates that a linear relationship between l/Pco 2 apP 
and v is obtained experimentally. It follows from hydrodynam- 
ics (Landau and Lifschitz, 1991) that the unstirred layer must be 
zero at zero viscosity. Therefore, the extrapolation of the linear 
regression line in Figure 15 to v = 0 gives the true membrane 
CO2 permeability Pco 2 true , which turns out to be 0.16 cm/s. The 
lowermost experimental data point was obtained in saline (arrow 
in Figure 15) and had a value of 0.12 cm/s. This means that the 
unstirred layer reduces Pco 2 by 25%, which can be considered a 
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minor effect. Inserting numbers into Equation (9), we can write: 
1/(0.12 cm/s) = 1/(0.16 cm/s) + 8/(2 • 10~ 5 cm 2 /s), 

and obtain an unstirred layer thickness 8 of 4.2 • 10~ 5 cm = 

0. 42 |xm in this case. This again is a very moderate unstirred layer 
thickness, which was obtained by optimized conventional stirring 
at maximal speed with a magnetic bar stabilized by being fixed 
in a Teflon ring. In conclusion, the Pco 2 values obtained by the 
mass spectrometric technique are affected by unstirred layers to a 
minor extent only. 

The empirical linear relationship of Figure 15 together with 
Equation (9) leads to the relation 

8/D C o 2 =C-v, (10) 

where C is a proportionality constant. This implies that 

8~v, (11) 

a prediction derived by Endeward and Gros (2009) from hydro- 
dynamic theory for this experiment, provided the magnetic stir- 
rer exerts a constant force in the cell suspension, independent 
of viscosity. It should be noted that this latter consideration 
merely provides a theoretical rationale for the empirical finding 
of Figure 15. 

A topic, although of great interest, has not been discussed here, 

1. e., the intracellular unstirred layer, which is omnipresent in all 
cellular transport measurements, but is not accessible experimen- 
tally. Endeward and Gros (2009) have used a computer model to 
estimate the size of this layer. For the case of CO2 diffusion in red 
cells they came up with a thickness of this layer of about V* to 1/3 
of the half-thickness of the erythrocytic disc. 

In view of the numerous vesicle experiments discussed in this 
article, it is of interest to finally add another insight from hydro- 
dynamic theory, which has been presented by Endeward and 
Gros (2009). When other conditions such as convection are held 
constant, the thickness of the extracellular unstirred layer is pro- 
portional to the square root of the diameter (or "characteristic 
length") d of the cell or vesicle studied 

8~Vd. (12) 

Using 8 = 0.42 \im for red cells with a diameter of 7.5 (xm, one 
predicts a 8 of 0.06 (xm for a vesicle diameter of 150 nm, i.e., 
the diffusion resistance of the unstirred layer toward CO2 will 
be seven times smaller in vesicles than in red cells. Thus, the 
estimated vesicular CO2 permeability will be corrupted by the 
unstirred layer to a much lesser extent. Applying an identical con- 
sideration on the basis of Equation (12) to MDCK cells (d = 
10 (Jim), one comes up with an estimate of 8 for MDCK cells of 
0.48 u,m, just slightly greater than for red cells. 
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